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Big bang simulation in superfluid ^He-B — Vor- 
tex nucleation in neutron-irradiated superflow 
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We report the observation of vortex formation 
upon the absorption of a thermal neutron in a ro- 
tating container of superfluid ■^He-B. The nuclear 
reaction n -|- ^He = p -|- fH -|- 0.76MeV heats a 
cigar shaped region of the superfluid into the nor- 
mal phase. The subsequent cooling of this region 
back through the superfluid transition results in 
the nucleation of quantized vortices. Depending 
on the superflow velocity, sufficiently large vor- 
tex rings grow under the influence of the Mag- 
nus force and escape into the container volume 
where they are detected individually with nuclear 
magnetic resonance. The larger the superflow ve- 
locity the smaller the rings which can expand. 
Thus it is possible to obtain information about 
the morphology of the initial defect network. We 
suggest that the nucleation of vortices during the 
rapid cool-down into the superfluid phase is simi- 
lar to the formation of defects during cosmological 
phase transitions in the early universe. 

An issue of current interest is the generation of topo- 
logical defects during symmetry-breaking phase transi- 
tions ||l| . In our experiment we study the density of vor- 
tex lines produced during a temperature quench of nor- 
mal '^He liquid into its superfluid B-phase. A sufficiently 
rapid transition results in a dense network of vortex Unes 
In a large enough bias field of superfluid flow some 
vortices will be stabiUzed and survive the subsequent re- 
laxation towards equilibrium while others contract and 
annihilate. Current theories of condensed matter physics 
generally provide the tools for understanding defect nu- 
cleation in slow phase transitions close to equilibrium. 
The present nucleation process, however, is far out of 
equilibrium and resembles the discussion of cosmologi- 
cal phase transitions in the early universe [^. Primor- 
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dial defect nucleation may explain the cosmological large 
scale structure and also in some scenarios the net baryon 
asymmetry of the present universe. To test the validity 
of such estimates it is important to find laboratory ex- 
periments in which defect formation in a rapid quench 
can be studied. 

Liquid '^He has several advantages over other systems 
which in the past have served as laboratory models for 
the cosmological defect-formation scenario, such as liq- 
uid crystals 0|, liquid ^He and superconductors [||. 
Many direct parallels and formal analogies connect su- 
perfluid ■'He theory with field theories which are used to 
describe the physical vacuum, gauge fields or fermionic 
elementary particles j7|. Superfluid ^He exhibits a vari- 
ety of phase transitions and topological defects which can 
be detected with NMR methods with single-defect sen- 
sitivity Of particular relevance to our rapid quench 
experiment is the fact that liquid "^He can be locally effi- 
ciently heated with thermal neutrons (Fig. |^). 

The neutrons are produced with a paraffin moderated 
Am-Be source of 10 mCi activity. They are then inci- 
dent upon the '^He sample container. At the minimum 
distance of 22 cm between source and sample, j/ « 20 
neutrons/min are absorbed by the '^He liquid. For ther- 
mal neutrons the mean free path is 0.1 mm in liquid '^He 
and thus all absorption reactions occur close to the walls 
of cylindrical container which has radius i? = 2.5 mm. 
Since the incident thermal neutron has low momentum, 
the 573 keV proton and 191 keV triton fly apart in oppo- 
site directions, producing 70 and 10 ^m long ionisation 
tracks, respectively. The subsequent charge recombina- 
tion yields a heated region of the normal liquid phase 
which for simplicity we here assume to be of a spheri- 
cally symmetric shape . 

The bubble of normal fluid cools by the diffusion of 
quasi-particle excitations out into the surrounding super- 
fluid with a diffusion constant D vpl where vp is their 
Fermi velocity and I the mean free path. The difference 
from the surrounding bulk temperature Tq as a function 
of the radial distance r from the centre of the bubble is 
given by 
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where Eq is the energy deposited by the neutron as heat 
and Cy the specific heat. The maximum value of the 
bubble radius Rb with fluid in the normal phase, T(r) > 
T., is 
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The bubble cools and shrinks away rapidly with the char- 
acteristic time Tq ~ Rb/^ ^ 10~^ s. The rapid super- 
fluid phase transition leads to the formation of a random 
network of vortices as was originally formulated by Kib- 
ble ||l| and later refined by Zurek ||]. 

The fate of nucleated vorticity depends on the bias 
field, the velocity Vg of the superflow outside the bubble 
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in the bulk ^He-B. The superflow is created by rotat- 
ing the container: the normal component is clamped to 
corotation with the container while the superfluid compo- 
nent remains stationary. This results in a relative super- 
fluid velocity Vg — at the wall of the container, when 
rotated at an angular velocity fl. If the radius of the 
vortex loop exceeds the value ro(ws) = (K/47rWi) Inro/^, 
where k = Trh/m^ is the circulation quantum and ^ = 
^0(1 — T/Tc)~^^^ is the superfluid coherence length, the 
loop expands. An expanding vortex ring eventually re- 
sults in a rectilinear vortex line which is pulled to the 
center of the container. If several vortex lines are formed 
they accumulate in a vortex cluster ||^ . 

The number of the vortex lines is monitored with 
NMR. In the insert of Fig. |two NMR absorption records 
are shown, starting from the moment when the neu- 
tron source is placed in position. The absorption events, 
which lead to nucleation, are visible as distinct steps. 
The step height gives the number of vortex lines nucle- 
ated per event, if the neutron source is sufficiently far 
from the sample and the absorption events are well sepa- 
rated in time. The number of nucleated vortex lines per 
unit time increases rapidly with increasing Vs (Fig. |^). 
The extrapolation to zero gives the threshold value Vcn 
plotted in Fig. ^ as a function of T for different pressures. 
This Vcn is smaller and has a different dependence on T 
than the critical velocity Vc at which a vortex is nucleated 
in the absence of neutrons js) . 

Let us estimate the critical velocity Vcn and the number 
of vortices Af{vs) created in one neutron absorption event 
The topological defects are nucleated during the 
superfluid transition to ^He-B when the order parameter 
begins to fall from the false ground-state into the true 
degenerate ground-states which form independently in 
regions which are not causally connected. The different 
regions grow in size and ultimately the order parameter 
fills all space but with a multitude of defects. These 
are laid down during the freeze out time ^/tqto , where 
To — ^o/vp characterizes the relaxation time of the order 
parameter in the low temperature limit. The value of 
the coherence length ^ at that point, Co('''q/t'o)^^'* [||, 
determines the initial distance between the defects, which 
in this experiment is of order 1 /im. 

The later evolution of the vortex network leads to 
a gradual increase of the intervortex distance 
caused almost entirely by the interaction of the defects 
with each other, rather than by fluctuations in the order 
parameter field as before. With increasing time loops 
with a line length I < ^ are smoothed out while larger 
loops straighten and reconnect such that the network ap- 
pears scale invariant. For a homogeneous system the 
number of loops n{l) per unit length and unit volume 
with line lengths I > ^ is given by 
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ical simulation gives the same distribution law but with 
slightly different C ^ 0.3. For the average straight-line 
dimension I? of a loop wc find V — PH^Y^^ with /3 « 1. 
The loop size distribution in terms of V is scale invari- 
ant, n{'D) dD « 2C ctD/V^, and does not depend on the 
distance between vortices. The evolution of the network 
leads to an increasing lower cutoff of the distribution, 
T^min = C(£, (t) , where a is of order unity. The upper cut- 
off is the diameter of the bubble 2?max = 2i?f,. 

When the average radius of curvature, determined by 
exceeds ro(ws) the vortices start to escape from the 
bubble. The number of vortices Af{vs), created per one 
neutron, is thus the number of loops with aro{vs) < 2? < 
2Rb within the bubble volume Vt: 



2Rb 



dV n{V) = ^ 



2Rb 



- 1 



(4) 



The critical velocity v^n is determined by the requirement 
U{vcn) = 0, which gives Vcn = ln(i?f,/^). Our nu- 

merical simulation suggests a ~ 1. Thus Vcn is inversely 
proportional to Rb and according to Eq. (0) has the tem- 
perature dependence Vcn oc (1 — Tq/TcY/ , in agreement 
with the measurements in Fig. ^. In terms of 
obtains the universal curve 
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where C ~ 0.4. For our case of finite volume of the 
bubble, where all vortices form closed loops, our numer- 



Eq. (H) reproduces the observed cubic velocity depen- 
dence and gives the correct order of magnitude estimate 
oiM{vs), as seen in Fig. |^. It also carries the same uni- 
versality feature as the measured results in Figs. Hand^: 
the ambient measuring conditions depend on what val- 
ues are chosen for the temperature T, pressure P, and 
magnetic field H, but in the results for N{vs) all this 
dependence is contained in the single parameter Vcn- 

The experiment demonstrates that defects are nucle- 
ated within bulk liquid '^He during a rapid quench to the 
superfluid state. An external bias flow, which exceeds 
a critical value, allows the defects to escape, to expand 
and to reach a stable state in the rotating container so 
that they can be detected one by one. With increasing 
superflow velocity smaller loops, which represent an ear- 
lier stage in the evolution of the network, are extracted. 
This allows us to probe the morphology of the defect net- 
work at different stages of its evolution. The measured 
number of extracted vortices as a function of velocity 
is consistent with an initial loop size distribution which 
characterizes the random phases of the order parameter 
in the quench. We expect to learn more about this distri- 
bution from measuring the velocity dependence of events 
which lead to a given number of extracted vortex lines. 
More sophisticated numerical simulation is expected to 
reveal how the network of large loops, which exceed crit- 
ical size, is modified by the bias superflow in the late 
stages of expansion. 
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FIG. 1. Principle of the experiment and of vortex nucle- 
ation during neutron irradiation: a) Superfluid ''He-B con- 
tained in a cylinder, which rotates at constant angular ve- 
locity n, is exposed to neutron radiation, b) In a neutron 
absorption event a proton and a triton are formed which fly 
apart in opposite directions. Their kinetic energies are used 
up in ionization reactions of ^He atoms, c) The ionized par- 
ticles recombine generating heat in the form of quasiparticle 
excitations which corresponds to a substantial part of the to- 
tal reaction energy. A cigar-shaped hot region of liquid is 
formed where the temperature rises above Tc. d) The hot 
bubble cools rapidly towards the bulk liquid temperature To- 
While passing through the superfluid transition at Tc a tan- 
gled network of vortex lines is formed. Vortex loops which 
exceed a critical size start to expand in the superflow created 
by the rotation of the container whereas smaller loops con- 
tract and annihilate, e) An expanding vortex ring ultimately 
intersects the walls of the container and gives rise to one rec- 
tilinear vortex line. The vortex lines accumulate in the center 
of the rotating container where they are counted with NMR. 



FIG. 2. Nucleation of vortices during neutron irradiation. 
Insert: NMR absorption as a function of running time t at 
high and low rotation velocity. The neutron source is turned 
on at t = 0. Each step corresponds to one neutron absorp- 
tion event. The height of the step gives the number of nu- 
cleated vortex lines and is denoted by the adjacent number. 
The high velocity trace is recorded with the neutron source 
at a 2.8 times larger distance from the '^He-B sample. The 
bulk liquid temperature is To — 0.96 Tc. Main frame: The 
rate N at which vortex lines are created during neutron ir- 
radiation, plotted as a function of the normalized superflow 
velocity Vs/vcn- Below the critical threshold at Vcn the rate 
vanishes while above the rate follows the fitted cubic depen- 
dence N = j[{vs/vcnf - 1], with 7 = (1.37 ± 0.03) vortex 
lines/min, shown by the solid line. The rate of vortex nucle- 
ation is the number of vortices per one neutron M{vs) mul- 
tiplied by the neutron flux v. N = vj\f{va). The theoretical 
estimate for Nivs) from Eq. (^ with v ^ 2Q neutrons/min 
gives 7 = v-kC 1% ~ 2, which is in order of magnitude agree- 
ment with the experiment. The filled data points are mea- 
sured at a bulk liquid temperature To = 0.96 Tc and the open 
points at 0.90 Tc- Both sets of data fall on the same cu- 
bic dependence although Vcn differs by 50 % (Fig. ^). The 
cubic dependence on the superflow velocity is ultimately in- 
terrupted by the spontaneous critical velocity Vc in Fig. ^, 
which is the maximum possible superflow velocity. The liquid 
pressure P = 18.0 bar and the magnetic field H = 11.8 mT 
are the same as in the insert. 

FIG. 3. Critical superflow velocities as a function of the 
normalized bulk liquid temperature 1 — To/Tc and pres- 
sure P of the liquid ^He-B sample in a magnetic field of 
H — 11.8 mT. Solid lines: Critical value of superfiow ve- 
locity Vcn oc (1 — T /TcY^^ at which vortex creation starts in 
the presence of the neutron source. Dashed lines: Critical 
threshold for superfiow Vc oc {1 — T/TcY^"^ in the absence of 
the neutron source ^]. 

FIG. 4. Scale invariance of the vortex nucleation law during 
neutron irradiation in different ambient conditions: in Eq. (0) 
the dependence on temperature, pressure, and magnetic field 
is contained in the critical velocity Vcn- In this plot the nu- 
cleation rate N is shown as a function of uj. The intercept 
with the horizontal Eixis gives Vcn for each line while the com- 
mon intercept with the vertical axis, 7 = (1.46 ± 0.12) vortex 
lines/min, defines the numerical prefactor in Eq. The ex- 
perimental variables P, T and H are listed in the legend for 
each set of data points. 
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